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ABSTRACT 



The collective excitation of hybrid atomic-molecular condensate are studied using variational 
method. The dipole response of the system is studied in detail. We found that the out-of-phase 
dipole oscillation frequency increasing slowly with the detuning parameter and is only few times 
large than the trap frequency making the experimental observation of this mode of vibration 
feasible. 
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The recent realization of Bose-Einstein condensation (BEC) in a dilute gas of alkali atoms 
has open a new opportunity for the theoretical and experimental investigation in quantum 
degenerate many-body systems [0, |[. In contrast to superfiuid helium, these weakly interacting 
gases are much more amenable to the theoretical prediction and quantative analysis. After the 
observations of such condensates, a new generation of more complex experiments has been 
done involving the production of multiple species BEC. Such mixture may consist of different 
hypefine states of same atom or different spin state of same stomic species. In the former case 
the composition of mixture remains fixed and in the later case the internal conversion leads 
to important new phenomena. It has been proposed that the interation that bring a binary 
atom system to an intermediate state molecule in the Feshbach resonance, recently observed 
in BEC0, ^, create in the dilute atomic Bose-Eintein condensate (BEC) a second molecular 
condensate component 0. The communication between the atomic and molecular condensates 
ensues through a Josephson-like oscillation of atom pair. It was found that the ground state of 
this hybrid condensate is a dilute BEC with liquid-like properties of a self-determined density. 
These findings were obtained using a homogeneous gas approximation (no confining potential). 
In this work we investigate the collective excitations of the hybrid system of atoms and molecules 
subjected to an appropriate trap. In particular we derive an appproximate expression for the 
dipole frequency of the out-of-phase oscillations of the atomic and molecular condensate centers 
of mass and investigate its dependence on the detuning and Feshbach resonance parameters. 
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The mean-field hamiltonian density describing the trapped system is 
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where V^a{m) denote the macroscopic wavef unction of the atom (molecule) condensate, e is the 
detuning parameter, a is the Feshbach resonance parameter introduced in Ref. 0, and we adopt 
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here a simple form for the atomic (molecular) confining potentials, Va = |^(A"x^ + A^y^ + A"^;^), 

2 

Vm = 2{-2m) (-^"-^^ + '^^v'^ + -^r^^)- usual, the Gross-Pitaevskii equations of motion for 
this model can be derived from the action T = J dt d?r^ i^ai^a — i^ai'a + — i'a'^a — 

J dt E[tpa, 'ipm] , where E[tpa, ^m] = Jd^rH . 

In order to investigate the general properties of the hybrid system described by H, Eq.(l), 
we employ the time- dependent variational Gaussian ansatz |H| 
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where qj,pj,rj,7ij {Qj, Pj, Rj,Il'^) are our variational parameters for the atomic (molecular) 
condensate wavefunctions and we have introduced the notation q = (qxQyQz)^^^ and Q = 
{QxQyQzY^^ ■ Thus, i/jai'ipm) is a Gaussiau centered at (Rj) with a width qj (Qj) and conjugate 
momenta 7r| (11^) and pj (Pj) respectively. In addition, we shall use as variational parameters 
the numbers of the atoms, / d^r\'il)*\'^ = Na, and molecules, / d^rl'ip^l'^ = Nm, subject to the 
constraint A^^ + 2A^m = (fixed). 



The action in this approximaion is therefore 

(4) 

The total energy has the general form E = Ea + Em + Ei + -E^^, where Ea{Em) are the 
usual Gaussian Gross-Pitaevskii energy for a single atomic (molecular) condensates 0, E'a = 

Nai:,=.,y,-Ah^] + M + £^7' + f + f<^|] + - n, and etc for 

Era)- The other two terms are, respectively, the elastic and inelastic coUisional contribution 
between the condensates, Ei = ^^^{nj=x, y, z exp[- ^^I^ ]{q^ + Q])-^/^}, E^^ = eN^ + 
2V2^3/4q3Q3 {^^^^ nj=x, y, z ^^P[jt' ~ '^ilS' ^''^ + c.c.}, where we have used the notations: Qj = 
B* + 2pj, hj = -2B*Dj - Ap^r^ + iU^ - 2nj, cj = B*D] + 2(3j{r^f - iUjU'^ - 27rjT|, being 
ttj = B* + 2(3, 2B = Q^2j — PjQ^^, 2(3 = q~2j — PjqJ^. The variational equations are then 
obtained from (^F = leading a set of first order differential equation of the form 

Xk = , (5) 

k * 

where S is an appropriate matrix. These classical equations of motion give a simple picture 
for the time evolution of the atomic (molecular) condensate wavefunction parametrized by its 
widths and centers of mass and the dynamics of inter-condensate tunneling driven by the FR 
interaction. The ground state energy, E = J d^r Ti , on the other hand is determined by the 
solution of Eq.(5) with Xk = 0. The small amplitude oscillation dynamics is described by an 
expansion around the equilibrium solution. Thus, if X^^^ is the static equilibrium solution and 
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X^^^ is the perturbation, we have 
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where Aj,- 



and K is an appropriate 26 x 26 inertia matrix. The set of 26 cou- 
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pled equations is diagonahzed for the frequencies, which come out in pair of opposite parities. 
The thirteen physical frequencies can be grouped into three sets: one corresponds to the mo- 
tion of the centers of mass of atomic/molecular gases; one defines the spatial frontiers of the 
atomic/molecular condensate and one describes the Josephson tunnelling effects between the 
two condensates. The frequencies associated with the Josephson tunneling were found in Ref . 
(no trap) to he u{Xa = ^,n) = ±nA[4Xa + Xay/l - X„ + y^]- 

An important new phenomenon in the hybrid condensate physics appears when the system 
is subjected to the confining trap namely the in-phase and out-of-phase dipole oscillations of 
the two centers of mass, what results four excitation modes when cylindrically symmetric trap 
is used. The dipole oscillation of the center of mass of an atomic BEC is entirely determined 
by the frequency of the trap and does not depend on the interatomic interaction 0. If, on the 
other hand, we allow atoms and molecules, with their own confining traps and with Josephson 
type oscillations, then presumably the dipole mode of oscillation of the respective center of 
masses will depend on the Feshbach resonance interaction parameter, a, and on the detuning 
parameter e. If so, then there is an interesting phenomenon to be explored experimentally. 
An analogous situation is discussed in a homogeneous hybrid condensate (no trap). In that 
case, the low k-modes are, as experted for superfluids, longitudinal waves, i.e. the oscillating 
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particle current is in the direction of propagation (k). One mode, in which both the atomic and 
molecular particle current move in the same direction, i.e. a Goldstone mode and its excitation 
energy vanishes as k — > 0. The other mode has atomic and molecular currents moving in 
opposite direction like a plasmon mode and it has a gap: it goes to a finite value which depends 
on the interaction parameters as k — > 0. In our case we would expect two similar modes: one 
gives an energy u>trap if atoms and molecules experience the same trap potential as in the case 
of the single condensate and another mode that has higher energy, which depends mostly on 
the interaction parameters. In the following we quantify these observations using our model. 

For simplicity we consider first the case of isotropic trap, A^, = A^^ = = 1, for both atoms 
and molecules. In this case, the following trivial equilibrium solution is found 

r, = TTe = i?c = n, = . (7) 

In addition, it is easy to show that q^^qx — where X*^ are generic coordinates of the 
subspaces corresponding the centers of mass and X'^ represent the other degrees of freedom of 
the variational space. Thus, the matrix A reduces to the form ^ ^ ^ ^ , where Cik = qx^qx'^ 
and Wik — QxwQxw ■ As a consequence, Xc decouple from Xyj at X° and the dynamics of the 
centers of mass system is described by two coupled oscillator equations. Eliminating tt^^ and 
n(o) we find 

Nat = -[Naiol+ Auj']r^ + Auj' (8) 
2N^R^ = -[2NmCol + Au;^]R^ + Au;^r^ (9) 

where the frequency shift Aa; is produced by the interaction between the atomic and molecular 



condensates and thus it contains the a— and e— dependence. This is related to the curvature 

of i?* at by the following simple expression 

A 2 ^ (d'^E\ (N N sfd^E\fd^E\ 

= "^"^(^) + ;^%) + (]^ + 2i^)(^)(^) 

The two coupled equations (8) and (9) can be uncoupled in favor of the normal modes, r] = 

NaTc + 2NmRc and ^ = rc — Rc, whose frequencies are, respectively. 

The mode f2_ is the in-phase solution where the centers of mass of atomic and molecular gases 
move in the same direction and is independent of the interactions 0. On the other hand, the 
out-of-phase dipole oscillation solution depends strongly on the composition of the system, 
which is self-determined by the equilibrium condition Xq. From the previous discussion we 
have ^ ~ NaN^ and ~ NaN^\ Thus, when « N, one has fi+ N^'^Z^ This 

result is illustrated in Fig.l, where f2+ increases dramatically for solutions with small number 
of molecules. 

We turn next to the quadrupole response of the system. In Fig. 2 we show the three corre- 
sponding frequencies The lowest frequency describes the oscillations q, Q and Na in the same 
direction. The other two modes describe the out-of-phase solution of q against Q, which leads 
to migration of atomic condensate to molecular condensate (dash-dotted plot) and vice-verse 
(dotted curves ). This can also be represented by the current ja = V Pa [jm = Vpm) obeying 
continuity equation V ■ {ja + 2jm) = 0. These modes, again, are strongly dependent on the 
particular equilibrium point of the system. 



In conclusion, we have discussed in this paper a variational calculation of the small amplitude 
dynamics of a hybridatom/molecule Bose Einstein condensate. The out of phase dipole mode of 
oscillation of center of mass of the two condensed gases is considered in detail. For reasonable 
values of the physical parameters (the Feshbach resonance strength parameter, densities of 
gases, etc), we have found that the dipole frequency depends approximatelly quadratically on 
the detuning parametes, e, at small value of the latter attaining an almost linear dependence 
for higher values of e. This frequency is less than an order of magnitude large than the trap 
frequencyfor values of e of several KHz. 
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Figure Captions 

Figure 1. The out-of-phase dipole frequency as function of detuning parameter, e, for the 
hybrid condensate in unit of trap frequency. The parameters are: Xa — Xm — ^ — 100m/x^s~l, 
N = 10^ K = 27r * 50Hz, a = 200m//3/2g g^^ ^^^^ details. 

Figure 2. The three eigenfrequencies for the q, Q and degrees of freedom. Same parameters 
are used as in Fig.l. Lower curve, uj{q + Q + No); middle curve, — Q + No); upper curve, 
u!{q — Q — Na). See text for details. 
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Figure 1 
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Figure 2 
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